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We demonstrate that the condensed matter quantum systems encompassing two reservoirs con- 
nected by a junction permit a natural definition of flows of conserved measures, Renyi entropies. 
Such flows are similar to the flows of physical conserved quantities such as charge and energy. We 
develop a perturbation technique that permits efficient computation of Renyi entropy flows and an- 
alyze second- and fourth order contributions. Second-order approximation was shown to correspond 
directly to the transition events in the system and thereby to posess a set of "intuitive" features. The 
analysis of fourth-order corrections reveals a more complicated picture: the "intuitive" relations do 
not hold anymore, and the corrections exhibit divergencies in low- temperature limit manifesting an 
intriguing non-analytical dependence of the flows on coupling strength in the limit of weak couplings 
and vanishing temperatures. 
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I. INTRODUCTION 

Modern condensed matter theory borrows concepts 
from quantum information theory, n particular, entan- 
glement. There is an explosive growth of applications of 
these conceps [l|]. Although the measures of entangle- 
ment are not linear in density matrix and as such cannot 
be directly related to quantum observables, these mea- 
sures proved to provide insight into complex structure of 
many-body wave functions and in many cases can be ef- 
ficiently computed numerically. [2^ Most applications are 
restricted to an equilibrium situation where the object of 
interest is the wave function of the ground state of a sys- 
tem 1^ >. The system is cut in to two subsystems A, B 
(typically, in real space). Partial trace over a subsystem 
gives a reduced density matrix BS^^ = Tr^{|^ >< ^|} 
in Hilbert space of another subsystem. Roughly, the von 
Neumann entropy of R can be ascribed to the bound- 
ary between the subsystems. There is a number of illu- 
minating theoretical results jl|, Q related to the von 
Neumann entropy of the boundary and its scaling with 
the boundary area. Not only entropy is under scrunity: 
the attention is paid to and the information is gathered 
from the whole spectum of the reduced density matrix 
R. [2I, [5I, Q The information about the spectrum can be 
compactified to a set of quantities that have been in- 
troduced in the context of information theory by Alfred 
Renyi and are called Renyi entropies. We define them 
as follows: Isl 



Sm = Tt 



Logs of Sm are obviously extensive quantities: if R can 
be presented as a direct product over a number of sub- 
systems, their contrubutions to In Sm add. The von Neu- 
mann entropy is obtained by taking a formal limit 

^ OSm dlnSM 

S = - lim — — - = - lim — 77T1—- 

M^i dM M^i dM 

Several works address time evolution of entropies in the 



course of transition processes like a quantum quench. [9|. 
There are interesting results concerning von Neumann 
entropy production in stationary non-equilibrium sys- 
tems, in particular, in quantum point contacts [Io|. A 
simplest quantum point contact is a single-channel con- 
ductor with transmission coefficient Tq, that connects two 
electronic reservoirs. The entropy flow for zero tempera- 
ture in the reservoirs was shown to be ^Ol 



dS 
~dt 



dt 



(TolnTo + (l-To)ln(l-To)). (1) 



Here, {dNg,tt/dt) is the number of electrons per unit 
time that attempt to transfer the contact, {dNg,tt/dt) = 
eV I (2i\K) The energy flow to either reservoir is given by 
dEldt = {dN^tt/dt)ToeV/2. The approach of |10| was 
rather heuristic relying on the representation of the elec- 
tron many-body state in terms of a sequence of individual 
single-electron scattering events [11], that is instrumen- 
tal in the field of full counting statistics. Levitov and 
Klich have rederived these results with more microscopic 
approach, extented those to the case of a general time- 
dependent scatterer and discovered a remarkable corre- 
spondence between the entropy flow and full counting 
statistics of electron transfers. This has been further 
elaborated in [l3|. 

There are two interesting peculiarities in the relation 
([1]) that have not received a proper discussion so far. 
Firstly, non-analiticity of the expression ([1]) at small Tq 
should indicate an intriguing divergence of perturbation 
series in Tq that measures coupling strength between the 
reservoirs. Secondly, second law of thermodynamics re- 
lates the entropy and heat increments for a system in 
thermal equilibrium at temperature T. From this one 
would conjecture the relation between entropy and en- 
ergy flows, 

dS _ 1 dE 

'dt ~ k^~dt' ^ ^ 

that would seem to have a textbook status. However, the 
result ([1]) is not compatible with this conjecture. The ini- 
tial motivation for the research presented in this Article 
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was to understand these pecuharities in a general cal- 
culable framework rather than for a restrictive case of 
quantum point contact. 



The main result of the Article is that the flows of Renyi 
entropies (Re-flows in short) are well-deflned in standard 
condensed matter setups comprising two reservoirs con- 
nected by a junction. Importantly, they permit a detailed 
evaluation in the framework of a quantum perturbation 
theory developed here. This is in contrast to the entropy 
flow: albeit the latter can be readibly obtained by an- 
alytical continuation from integer to continous M. We 
study second and fourth order of perturbation theory. 
The second-order results were shown to be very special: 
they satisfy the "textbook" relation ([2]) as well as some 
extra "intuitive" relations outlined in the text. This is 
because the second-order results can be put into direct 
correspondence with the transition events taking place in 
the system. Neither this correspondence nor "intuitive" 
relations hold for forth-order corrections. These cor- 
rections exhibit non-analytical low-temperature behav- 
ior manifesting non-analytical dependence of the flows 
on coupling strength in the limit of weak couplings and 
vanishing temperatures. 



An obvious point of critisism is that the Re-flows are 
"unphysical" since they are non-linear in density matrix 
and therefore cannot be readily asscociated with (the 
flows of) measurable quantities. The author shares the 
opinion and idea of Levitov and Klich [12] that the flows 
of quantum infromation quantities may be related to the 
statistics of flows of quantum observables by universal 
relations whose general form is yet to discover. If this is 
true, a measurement of such statistics would in fact con- 
stitute a measurement of an "unphysical" flow. In any 
case, the Re-flows characterize information exchange be- 
tween the reservoirs in a rather detailed way. This makes 
their evaluation useful. 



The structure of the paper is as follows. In Section 
[m we detail and illustrate the deflnition of Renyi en- 
tropies. In Section lllli we show that Re-flows can be 
associated with conserving currents, these currents be- 
ing rather similar to electric or energy current through a 
junction in quantum transport setups. Next SectionllVlis 
devoted to the description of a generalized Keldysh tech- 
nique that enables perturbative calculation of Re-flows in 
an arbitrary setup. We evaluate the flows in the secon- 
order approximation in Section |Vl The second-order ap- 
proximation is characterized by a set of speciflc relations, 
those outlined in Sections I VI) I VIII In Section IVlIII we de- 
rive the fourth-order corrections. We demonstrate that 
these corrections diverge in the limit of vanishing temper- 
ature (Section lIXp indicating non-analytical dependence 
of the flows in the limit of small couplings. We conclude 
in Section 1X1 



II. CONSERVED MEASURES 

Let us discuss the deflnition and straightforward prop- 
erties of Renyi entropies. We start with an isolated flnite 
quantum system that is characterized by density matrix 
R. We define the Renyi entropies as traces of integer 
powers of 

5M = Tr{^^} (3) 

This definition can be easily extended to non-integer M, 

Sm =Tr{exp (mIu^)} 

this identifies Sm as Laplace transform of the spectral 
function of the operator In R. This identification permits 
to relate the von Neumann entropy and Renyi entropies. 



S = -TT{R\nR} 



hm Sm 



lim In^M. (4) 



If the system is in thermodynamic equilibrium at 
temperature T, the Renyi entropies are readily ex- 
pressed in terms of the temperature-dependent free en- 
ergy F(T) (e.g. [14]), 



In^M 



M 



{F{T/M) - F{T)) . 



(5) 



The Renyi entropies can be easily computed for simple 
non-equilibrium quantum systems as well. The example 
we will use in this article concerns free fermions in a 
system of single-particle levels labelled by k. We ascribe 
the levels arbitraty filling factors fk. The density matrix 
of the system reads 

k 

fik = a^a/e being number operator, ak being annihilation 
operator in the level /c, fk = 1 — //e- From this, the Renyi 
entropies read 



In 5m = In (/f + 



(6) 



For a degenerate Fermi gas in thermal equilibtrium 
with constant density of states 5^^ near the Fermi energy, 
this reduces to 



IuSm 



QSs 



77 
M 



(7) 



The quantum evolution of the system is governed by 
a Hamiltonian H. The same Hamitonian determines the 
evolution of the density matrix. 



-ih 



dR 
dt 



[H,R]. 
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Since the density matrices in different momets of time 
are related by unitary transform, the trace of any power 
of R does not depend on time. 

We just proved that Renyi entropies provide a set 
of conserved measures for density matrix of an isolated 
quantum system, 

—Sm = 0, 
at 

in similarity to the conserved physical quantities like en- 
ergy or change. The difference is that Sm are not linear 
in density matrix and therefore cannot be immediately 
associated with any quantum observables. 

While these definitions and properties may seem 
straightforward, there are some caveats to discuss. They 
are related to the traditional classical definition of en- 
tropy and its use. To illustrate, let us start with the 
non-equilibrium Fermi system in a pure state where fk 
are either 1 or 0. Obviously, Sm = 1 for this pure state. 
Let us take H that involves weak interaction between 
the fermions. Common knowledge suggests that the in- 
teraction causes thermalization of the fermion distribu- 
tion function. After some time, the filling factors fk will 
correspond to a Fermi distribution at some effective tem- 
perature T* that depends on the initial state. It may 
seem that the entropies are now given by Eq. [71 and thus 
differ from 1. However, they must conserve! Indeed, af- 
ter the Hamiltonian evolution the system is still in a pure 
state. 

The solution to this apparent paradox is that Eq. [6]can 
not be used for an interacting Fermion system whatever 
small the interaction is. We could have used an alterna- 
tive definition of the eutopics that is close to the original 
classical definition of Boltzman. Namely, we could reduce 
the true density matrix to 2 x 2 matrices pk by taking 
partial trace over all fermion Fock states not involving 
the level /c, take the diagonal elements of these matri- 
ces and substitute them to [6l Such alternative definition 
might even seem more useful, for instance, for description 
of thermalization in an electron gas. However, this defini- 
tion would necessary involve entirely subjective choice of 
a basis (that of non-interacting fermion Fock states) and 
would not give rise to conveniently conserved measures. 
On these grounds, we stick to "quantum" definition (|3|). 

Let us now consider a bipartition of the Hilbert space 
A ^ B corresponding to two systems A and B. We can 
now define two sets of Renyi entropies, 

5if)=Tr.{(i?(-))"};5if)=IV.{(i?(-))"}; 

(8) 

where the reduced density matrices in two subspaces are 
defined with the aid of the partial traces in these sub- 
spaces, 

R^^^ = TiBiR}; R^^^ = TrA{^}. (9) 

If the quantum evolutions of the systems are completely 
independent, 

H — Ha + Hb, 



Ha,b being operators involving the corresponding sub- 
spaces only, both sets provide the conserved measures. 

It is interesting to note that the sets of Renyi entropies 
are not the only conserved measures characteristic for a 
bipartition. Any polynomial in density matrix that is 
invariant with respect to the group Ua^Ub of unitary 
transforms in two subspaces, would provide such a mea- 
sure. To give a minimal example, let us label the states 
in A[B) with Latin (Greek) indeces. The quantity 

K = ^ Raa,b-fRb^,caRc-f,a(3 (10) 
a,6,c;Q;,/3,7 

is a conserved measure that can be reduced neither to the 
Renyi entropies of the systems nor to the Renyi entropy 
of the whole system. The characterization of all such 
measures forms an interesting research task beyond the 
scope of this Article. 

III. FLOWS 

Let us complicate the bipartition situation by including 
a Hamiltonian coupling between the systems A and 5, 

H = Ha^Hb^Hab. (11) 

Hab being an operator that involves degrees of freedom 
in both subspaces. 

We will assume that the systems A and B are in- 
fitely large and are characterized by continuous exci- 
tation spectrum given by Ha,b respectively. However, 
Hab couples a relatively small number of degrees of free- 
dom in both systems. It is convenient to assume that 
the systems A, 5are in thermal equilibrium at different 
temperatures Ta^Tb and/or chemical potentials fiA, Pb 
and thus play a role of reservoirs. The coupling is such 
that it does not lead to equilibration of Ta^Tb {i~ia-,I^b)- 
Rather, there is a constant energy /particle flow between 
the systems that does not depend on actual volume of A 
and 5, nor on their detailed properties. 

Examples of such arrangement are common quantum 
transport setups. [l5[ An exemplary setup (Fig. [1]) con- 
sists of two (infitite) current-carrying leads A and B kept 
at different chemical potentials. These leads are con- 
nected by a junction. Simplest coupling that we will use 
in our examples is the tunneling Hamiltonian 

Hab = ^^A,/CB i^l^dk^ +h.c) (12) 

that describes particle transfer between the leads. 

It is important in quantum transport setups that the 
value to measure experimentally and compute theoret- 
ically is a flow of a conserved quantity, electric charge 
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FIG. 1: Two leads (reservoirs) A and B in a typical quan- 
tum transport setup correspond to the bipartition A ^ B. 
The junction connecting the two is associated with coupling 
Hamiltonian H^^^^ . In similarity with the flows of charge and 
energy, one can define the flows of Renyi entropies (Re-flows) 
in such setups. 

being a simplest and most usefull example. In principle, 
one defines a flow of charge through a cross-section of the 
setup, this cross-section may also define the bipartition. 
Owing to conservation, this flow does not depend on the 
cross-section. Therefore, the exact bipartition is also not 
important. Owing to this, the current is said to depend 
on the properties of the junction rather than on details 
of the reservoirs. This enables the experimental investi- 
gation of nanostructures: a signal in a measurement of a 
non-conserving quantity would be most probably domi- 
nated by massive leads rather than by a small nanostruc- 
ture, while the flow of conserved charge is determined 
by the nanostructure forming the junction between the 
leads. The same pertains energy flow. Energy and charge 
are not the only conserved quatities: approximately con- 
served spin currents [l6| provide another practical exam- 
ple. Artificially constructed conserved quantities appear 
useful in description of quantum coherence in nanostruc- 
tures, being mathematical basis of so-called circuit theory 
of quantum transport. [l5| 

From analogy with the flows of conserved quantities 
in quantum transport setups, we conjecture that there 
are finite flows of conserved measures, Renyi entropies 
(Re-flows), to each subsystem A and 5, 

m Dj^ — j-j^ , 

Owing to conservation of Reneyi entropy in each system, 
the Re-flows would not depend on exact bipartition of the 
system and are determined by properties of the junction 
that is in principle described by Hab- 

There is an important difference. For physical quanti- 
ties the conservation holds in the whole system as well as 
in each subsystems. For instance, elecrtical currents to 
each lead must satisfy Ia-\-Ib = 0.[17] As far as Renyi en- 
tropies are concerned, there is no exact conservation law 
for a sum In^*^^ ^In^*^^ at finite Hab^ although these 
quatities are extensive. There is a conservation law for 
the total Renyi entropy \nS^^^^\ However, the latter 
at finite Hab is the sum In 6*^^ + In 6*^^ only approx- 



imately, up to the terms proportional to the volume of 
the system. 

Therefore, in general 

Let us compute these Re-flows. We will always restict 
ourselves to the Re-flow to the system A: the Re-flow to 
the system B is obtained by permutation of A and B. 
For brevity, we will skip the index A in s\^^ and 
where this does not lead to confusion. 



IV. PERTURBATION TECHNIQUE 

We will use a perturbation technique in Hab and keep 
the calculation as general as possible. We assume "adia- 
batic switching" of the perturbation[18]: far in the past 
the coupling is absent, and the density matrix is a direct 
product over subspaces, 

^(-oo) = Ra{-oo) (g) Rb{-oo) 
Ra{-oo) ^Pa\a >< a\; Rb{-oo) ^Pa\a >< a\. 

a a 

(As above, we label the states in subspaces A (B) with 
Latin (Greek) indexes. ) The coupling slowly grows 
achieving actual values at time t. The time evolution 
is given by (from now on, h = 1) 

R{t) = Texp (^i J drHABir)^ R{-oo) x 

xtexp(^-z J drHABir)^ 

Hab(j) is taken here in interaction representation, 
Texp(Texp) denote time(anti)ordering in the evolution 
exponents. Expanding this in HABij) gives perturba- 
tion series most conveniently presented as diagramms in- 
volving the Keldysh contour (Figl2]). The operators in 
perturbation series are ordered along the contour. Two 
parts of the contour correspond to time evolution of bra's 
and ket's in the density matrix. The crosses represent 
the (time-dependent) perturbation HABif) at a certain 
time moment. The intergatation over time moments of 
all perturbations is implied. There is a state index asso- 
ciated with each piece of the contour. Since R{—oo) is 
diagonal, this index does not change when passing this 
element. The index changes if a non-diagonal matrix ele- 
ment of the perturbation is considered. Summation over 
indices is implied. 

The crucial observation is that this scheme can be 
straightforwardly generalized to any integer number M 
of density matrices. These matrices undergo independent 
unitary evolution in time interval (— oo, t). It is construc- 
tive to think of a set of M "parrallel worlds" and draw 
the diagramms for perturbation series using M parallel 
bra- and ket-contours. To compute 6'm(^) with this set, 



5 




FIG. 2: Perturbation theory for a single density matrix on 
the Keldysh contour. 




T 



FIG. 3: A diagramm of perturbation theory for S\^^ for M — 
3. It involves three parallel worlds. Reconnect ion of Keldysh 
contours for subspaces A (black) and B (white) accounts for 
partial trace over B and matrix multiplication in A. 



we first need to 'split' the contours to account for pos- 
sibly different ordering of operators in subspaces Aand 
B (black and white curves in Fig. [3]). Then we need to 
reconnect the countours at r = All white contours are 
closed within each world, this corresponds to the partial 
trace over B for each density matrix involved. In contrast 
to this, the black contours are connected to form a single 
loop going through all the worlds, this corresponds to the 
matrix multiplication in the definition (|3]) of Renyi en- 
tropy. This conveniently represents the rules of operator 
ordering for any diagram of particular order in Hab- 

It is interesting to note that the reconnecting the con- 
tours in a different fashion gives rise to perturbation the- 
ory for other conserved measures, for instance, for K 
given by Eq.eqiother-measure. 

It is natural to require that the matrix elements of 
Hab are only non-diagonal, that is, i^^^f^ = if either 
a = b or a = p. In this case, the first non- vanishing 
contibutions to Re-flows will be of the second order in 
Hab- 



V. SECOND ORDER 

Let us compute the Re-flows in the second order in 
HaB' It is proficient to directy compute the time- 
derivative of Sm' For diagramms, this corresponds to 
placing one of the perturbations at r = t. The only way 
to satistfy the continuity of state index along the white 
contours is to place the second perturbation in the same 
world. Four contibuting diagramms are given in Fig. |Vl 
In fact, the same four diagramms arise in the derivation 
of Golden Rule transition rate and are familiar to anyone 
who made use of Keldysh perturbation theory for den- 
sity matrix. The specifics of Renyi entropies is reflected 
in extra factors p^~^ the diagramms acquire in compar- 
ison with the case of a single density matrix. Summing 
up the four diagramms yields 



a,a;b,(3 



(13) 



a,a;b,(3 



f dt'2Re U^t-t')(M^-E,-E,)\ 
The integral over time reduces to 

manifesting energy conservation between the initial state 
\aa > and final state 1 6/3 >. 

This suggests that we can rewrite the whole expression 
in terms of Golden Rule rates Taa,b^ of the transitions 
between the states \aa > and >, 

Taa,bp = 2T,\H['^ll\^5{Ea + E„ - E„ - Ep). (14) 
With this, the flow reads 

{Sm)J^M = M ^ Taa;bl3{PbP/3 - PaPa)P^~^ (15) 

a,a;b,(3 

We see that the flow vanishes if the systems are in ther- 
modynamic equilibrium at the same temperature. In- 
deed, in this case PbP^/PaPa = exp((£^5 Efs — Ea — 
Ea)/kBT) = 1. 

Since the transition rates T^aa 

,5/3 in Golden rule ap- 
proxiamtion are symmetric with respect to a permutation 
aa ^ b/3, we can regroup the terms to arrive at 

{Sm)Tm = Mj2^a^bPa{pf-' -P^-')- (16) 
a,b 



where 



a,f3 
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FIG. 4: Second order diagramms for time derivative of a 
Renyi entropy. The contibutions come only from perturba- 
tions in the same world, only this world is shown in 
each diagramm. For all diagramms, the perturbations are 
taken at time moments t and t' < t. The letters at the con- 
tours label the states involved. 

gives the total transition rate from the state |a > to the 
state \h > averaged over all possible configurations of sys- 
tem B. Let us use Eq. [16] to derive a simplified expression 
valid in zero-temperature limit. In this limit, the system 
A is initially in the ground state |0 >, so that po = 1 and 
Pa = for a ^ 0, Sm = 1- We obtain 

Tm = -MTo; (17) 

Fo being the total transition rate from the ground state 
to any other state. Remarkably, this involves no assump- 
tion concerning the system B: it can be very far from 
equilibrium. 

Eq. [16] is also a convenient starting point to derive the 
expression for the von Neumann entropy flow. Taking 
the limit M ^ 1, we obtain 

dS 

- ^ = ^ In (Pb/Pa) Ta^bPa- (18) 

a,b 

Let us assume thermal equilibrium of A. In this case, 
\n{pij/pa) = {Ea - Eb)/kBT). Summing up the energy 
changes £^5 — Ea in the course of individual transitions 



from a to 6, we prove that the energy flow to the system 
A equals 

^=Y.^a^b{E,-E,)p, 

a,b 

Comparing this with Eq. [181 we recover the relation (|2j) 
that appears to be universaly valid within the second- 
order perturbation theory. Remarkably, this invloves no 
assumption about the system B. 

One can easily imagine the situation where dE/dt re- 
mains constant at T ^ (see Eq. [1]). The entropy flow 
in this case diverges as 1/T. One can wonder how it is 
compatible with the fact that Re-flows approach finite 
limit at vanishing temperature. The point is a hidden 
non-analyticity of ([T7|) : the flow does not vanish in the 
limit M ^ 1 as it is implied by more general Eq. [16] 
Therefore, Eq. [T7] cannot be used to derive the flow of 
von Neumann entropy. 

VI. EXAMPLE: FERMI GAS 

Let us explicitly compute the Re-flows for tunneling 
between two non-equibrium Fermi gases. In this case, 
a (or b) is a number state defined by a set of occupa- 
tion numbers in all levels k {n/e}, Pa = Ylkifki^ — ^/e) + 
//c(^/c))- The perturbation Hab is given by Eq. [12] 
Let us concentrate on the transitions involving a level 
q. Since the probabilities factorize, we can forget about 
the Fock states involving all other levels since their con- 
tributions to right hand side of Eq. [16] cancels with their 
contribution to Sm- We can therefore consider only two 
states spanned by level subsituting into Eq. [16] their 
contribution to Sm = + f^- The probablilities of 
these states, |0 > {uq = 0) and |1 > {uq = 1) are fq 
and fq respectively. Substituting this into Eq. [16] and 
summing up over q yields 

^M = Mj:{Ttf,-T-f,) '^ , (19) 

q Jq ^ J q 

F+ (F~) being the rates of addition (extraction) of a 
fermion to (from) the level q. 

Let us outhne a heuristic way to obtiain the above rela- 
tion. This way does not involve any quantum mechanics: 
it is something L. Boltzmann would do if given the prob- 
lem. Let us treat In Sm and fq as continuous quantities 
related by Eq. [6] By virtue of the relation, the flux of 
In Sm would be expressed in terms of fluxes of fq as 




Since fq is the average number of particles in the level q, 
^ll — r+ f — r- f 

1. - ^ qJq ^ q Jq 
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Combining these two relations reproduces Eq. [191 We 
stress that we do not see any reason for this relation to 
hold in general, since both Sm and Uq change in discrete 
rather than continuous fashion in the course of a tun- 
neling transition. Indeed, the inspection of forth-order 
terms in perturbation theory shows that the heuristic 
way does not generally work. 

Let us assume that the filling factors in the system A 
depend on single-particle energy E only. Eq. reduces to 



M 



M 



I dE (r+ 



/-r-/ 



/ 



M-1 



/ 



M-1 



M 



(21) 



V^'~{E) being addition/extraction rates per energy in- 
terval. In case of vanishing temperature (in the system 
A), fq = 0{—E-\-fi). The Renyi entropy flow reduces to 

J='m=M J dE {T^{E)e{E - ii) 

-r-(E)e(-^ + /i))sgn(E-/i) = 

being numbers of added/extracted particles. This is 
in agreement with the general relation (pT|) since the only 
events that happen in our setup are indeed additions or 
extractions of single particles. As in general case, ana- 
lytical continuation to M ^ 1 at vanishing temperature 
does not work since the expression does not vanish in the 
limit M ^ 1. One can work at non-zero temperature, 
where it does vanish. Taking derivative of Eq. [21] at 
M ^ 1 yields the flow of von Neumann entropy 



dS_ 
dt 



j dE {r^{E)f{E) - r-{E)f{E)) In 



Substituting Fermi distribution, we recover once again 
the relation ([2]) between the entropy and energy flows. 



VII. MASTER EQUATION APPROACH 

In a single world, there is a simple way to resum the 
perturbation series and arrive to a master equation that 
contains only diagonal elements of density matrix (Fig. 
[SJ. For a diagramm, we split the time-line by perturba- 
tions into the blocks as shown in the Figure. The blocks 
come in two sorts: diagonal with the same state index in 
both contours and non-diagonal ones. Intergation over 
time duration of non-diagonal blocks gives Golden Rule 
transition rates. The integration over time duration of 
diagonal blocks is divergent, this indicates that the dia- 
gramms need to be resummed. The resummation gives 
the master equation invoving the diagonal elements of the 
density matrix only. In case of bipartition, this equation 
reads 



p{-oo) 




FIG. 5: Master equation in a single world is obtained from 
the resummation of perturbation series whereby the time- 
line is separated into diagonal and non-diagonal (grey-shaded) 
blocks. One can do this resummation in multiple worlds to 
compute the time evolution of Renyi entropies. 



Here Taa stands for the total transition rate from the 
state \aa >, that is the sum of the partial transition 
rates, F^q, = J^^^ ^aa^b/s- 

Inspired by the above exercises with second-order dia- 
grams, one can attempt the resummation in parallel worls 
in such a way that the perturbations in a non-diagonal 
block are situated in the same world. In this case, we have 
independent master equations governing the dynamics of 
density matrix in each world. 

Suppose we know the propagator of the master equa- 
tion that is, a linear kernel V expressing the evolution of 
the probabilities, 

Paa{t) = Vaa,h(3{t,t')pi,i3{t'). 

With this, we can easily express time evolution of Sm 
in time interval (0, t) 



M 



(23) 



a \ a 



and compute the Reneyi entropy flows. We will not an- 
alyze the general situation here since as all systems de- 
scribed by a master equation it belongs to the realm of 
classical physics anyway. Rather, we concentrate on a 
simple but constructive example. We will assume that 
initially the systems were in a pure product state. Let 
us denote this state |00 >. Let us further assume that 
the transitions from this state as well from all other state 
can happen to a great number of states. This number is 
so big that once the flrst transition has taken place the 
system will never get back to the original state. With 
this, the propability to remain in the state |00 > falls off 
exponentially. 

Poo = exp (-Foot) , 

Foo being total rate of all transitions from |00 >. Under 
same conditions, the sum in Sm is dominated by the 
probability to remain in the state |00 >. The SM{t) just 
equals the probability that during this time interval no 
transition event occurs in any of M parallel worlds. 



dt' 



^aaPa 



(22) 



SMit)=pZ=e^pi-Mroot). 



(24) 
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The flow is therefore Fm = — ^roo, in fuh correspon- 
dence with Eq. [13 

We note that the resummation of the diagramms we 
have done is not complete since we have assumed that 
each non-diagonal block comprises only two perturba- 
tions. The transition rates are thus approximated by 
their Golden Rule expressions. In most cases involving 
master equation in a single world, this is not important: 
the further orders of perturbation theory provide small 
corrections to the rates not changing the structure of the 
equation. It looks like this is not so in the case of parallel 
worlds. The analysis of fourth-order perturbation correc- 
tions presented below indicates the problems invloved. 

VIII. FOURTH ORDER 

Let us analyze the fourth-order diagramms for time 
derivative of Sm- As above, we assume that Hab does 
not contain diagonal elements. Since white contours are 
closed within each world, the four perturbations can ei- 
ther all come in the same world or in two pairs in two 
different worlds. If all four come in the same world, they 
describe a correction to one of the Golden Rule transition 
rates. This correction does not bring anything new and 
we disregard these diagramms in further consideration. 

A diagram involving two different worlds is given in 
Fig. IVIIII We see that in general the black contour en- 
tering a world with perturbations exits it with a different 
state index. For a particular case when these indices are 
the same, a = 6, the diagramm diverges upon integra- 
tion over time. This is not surprising since we expand 
SM{t) cx exp(J^M^)- The fourth-order expansion thus 
contains terms oc (J^^^)^t/2, J^^^^ being the second-order 
contibution to the rate that we have already calculated. 
Indeed, the diagram with a = bis proportional to (J^*^^^)^ 
and therefore does not contibute to fourth-order correc- 
tion to the flow. We thus concentrate on the case a ^ b. 
We call this diagram " quantum" since we will see that it 
does not permit an interpretation in terms of " classical" 
transition events. All expressions for T^dS/dt in this 
Section give fourth-order corrections to these quantities. 

There are 16 diagramms of this sort corresponding to 
the number of ways the pairs of H^^^^ in each world can 
be placed on bra and ket contours. Summing up all of 
them, we can present the fourth-order correction in the 
following form: 

^/^^ =^Y.\^ab?^{Ea- E,f-^ (25) 

, Pa Pb 

a,o 

A _ \^ tt{AB) rriAB) 

^ab — -^aa,c(3-^c(3,ba 
c,a,(3 

(tT {{Pa + Pb)Pa - '2pcPl3) S {Ea + E^ - E^ - Ep) 

-z 2^^^ y 

Ea ^ Ea — Ec — Ej3 J 

The structure of the matrix elements in the " amplitude" 




FIG. 6: A fourth-order "quantum" diagramm for Renyi en- 
tropy flows. The contibutions come from perturbations H^^^'^ 
in two different worlds, only these two worlds are shown. The 
letters on the contours label the states involved. 

Aab is the same as for an amplitude of the transition from 
the state \aa > to the state \ba >^ that is, without the 
change of the state of the subsystem B. Such transition 
would seem to involve a virtual state |c, /3 >. However, 
the rest of the expression for Aab does not support this 
interpretation: rather, probabilities enter in a form sug- 
gesting that the transition takes place between one of the 
states \aa >, |6a > and the state \c/3 >. Therefore, the 
expression can be associated with no "classical" transi- 
tion and corresponds to no actual transition rate. 

Let us assume that the probabilities in the system 
A depend only on energies of the corresponding states. 
Then it follows from Ea = £^6 that Pa = Pb- The term in 
Aab with the energy difference in the denominator van- 
ishes and the flow reduces to 

SmJ'm = (M - l)^^\Aab\^S{Ea - E,)p^-^; (26) 

a,b 

Aab - ZTT ^aa.c^^c^M 
{PaPa - PcP(3) S {Ea E^ - Ec - E^) . 

We notice that if both systems are in thermal equilib- 
rium, it follows from Ea-\-Ea = Ec^Efs that PaPa —PcPjS 
and the "amplitudes" Aab vanish. 

Let us compute the fourth-order contribution for non- 
equilibrium Fermi gas. We pick up two levels in sys- 
tem A, /ci^2, with Ek-^ = Ek^ = E. Their contibution 
to Sm reads Sm = {fA{E) + fA{E)Y. Let the state a 
correspond to |1 >k^ |0 and state b correspond to 
|0 >fei |1 >/c2* The probablilities of these states equal 
fA{E)fA{E). The "amplitude" Aab is contributed by 
all levels q in the system B with the energy Eq = E. 
There are two possible "virtual states" |c/3 >. First pos- 
sibflity is an empty level q and both levels ki^2 occu- 
pied (pa = fsiE), pp = ME), Pc = fl{E)). Sec- 
ond possibility is an occupied level q and empty levels 
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fci,2(pa = ME), VP = fsiE), Pc = flm. Owing to 
anticommutativity of fermionic creation/annihilation op- 
erators, these two contibutions to the "amplitude" come 
with opposite sign yielding 

Aab = MfB{E) - fA{E)) J2 tkuatk,,AEa - E). (27) 



To get the contibution of ki^2 to the Re-flow, we square 
Aah^ multiply it with p^"^ = {fAfA)^'^^ and divide with 
the contribution Sm = {Ia + IaY these levels to 
Renyi entropy. Next to it, we sum over all states in the 
energy interval {E^E -\- dE) and integrate over energy to 
obtain the total Re-flow. We obtain 



Tm = (1-M) j dE{fB-fA) 



where is deflned as 



2 ifAfAr-' 



EiE); (28) 



= StT^ ^ tk,qtk',q'tk',qtk,q' 

xS{E - Ek)6{E - E'^)5{E - Eq)5{E - E'^). 

IX. LOW-TEMPERATURE DIVERGENCIES 

We will now demonstrate that the "quantum" con- 
tibution derived manifests serious problems with term- 
by-term perturbation theory in the limit of vanishing 
temperature, indicating non-analytical dependence of the 
flows on coupling strength in the limit of weak couplings 
and vanishing temperatures. 

First let us note that the contibution seems to have 
an evident zero- temperature limit, namely zero, at least 
if the ground state of the system A is not degenerate. 
Indeed, delta-funcion in Eq. [25] cannot be satisfled for 
any state b ^ a. However, analythical continuation to 
non-integer M gives rise to problems. Let is illustrate 
this computing (|28|) assuming thermal distribution and 
the difference between chemical potentials in A and 5, 
ey, is much bigger than ksT. To this end, we can set 
= in Eq. Disgerarging the energy dependence of E 
near Fermi surface, we obtain 



M 



kuTE- 



l-M 



M 



1 



M Vsin^ 



27r 



(29) 



While the flow vanishes at vanishing temperature, it has 
a flnite limit = SksTE at M ^ 1 where it should 
vanish. 

This signals the divergence of fourth-order contribu- 
tion to the flow of von Neumann entropy. Indeed, from 



the general expression (|26|) we deduce that 

= y^\A,,\^5{Ea-Eh)-; (30) 



dS 
~dt 



Pa 



that is, the states with lesser probabilities Pa contibute 
most to the entropy flow! For Fermi gas example, the 
entropy flow reduces to 



dt J IaIa 



(31) 



For equal temperatures in the systems A and 5, this 
gives 



dS 



eV 



-— = AksT E sinh^ , 
dt 2k bT 



(32) 



The fourth-order contibution to the entropy flow thus 
diverges exponentally upon T ^ 0. 



X. CONCLUSIONS 

In this paper, we have shown that the flows of Renyi 
entropies are well-deflned for common setups encompass- 
ing a junction between two reservoirs A and B. They are 
similar to flows of physical conserved quanities like en- 
ergy and charge. They vanish in thermal equilibrium 
and can be induced by the difference of chemical poten- 
tials/ temper at ues in the reservoirs. The flows charac- 
terize (quantum) information transfer between the reser- 
voirs and generalize the flow of von Neumann entropy. 

We develop a general quantum perturbation technique 
to compute these flows, this involves a set of parallel 
worlds. We consider second and forth order terms. In 
second-order approximation, the flows can be related to 
Golden Rule transition events. As such, they satisfy a 
set of "intuitive" relations expressed by Eqs. [21 [231 [201 

Our analysis of fourth-order corrections reveals in- 
teresting physics: the "intuitive" relations appear to 
be wrong, and corrections exhibit divergences in the 
limit of low temperatures manifesting an intriguing non- 
analytical dependence of the flows on coupling strength in 
the Hmit of weak couplings and vanishing temperatures. 
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